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Neutron kinetics 


This application note deals with a simple implementation of neutron kinetics for 
nuclear reactors fuelled with 7?°U or 73°Pu respectively. It will be part of a little series 
of application notes culminating in the simulation of a nuclear power plant on an analog 
computer. 

The neutron kinetics simulation is a central part of every nuclear reactor simulation.* 
Following [FENECH et al. 1973, p. 128] neutron kinetics can be described by 
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Here, 7 is the neutron density, s represents an external neutron source (which is 
often used to “start” a nuclear reactor), /* represents the effective neutron lifetime, 
which depends on the reactor design and geometry and can range between values as big 
as 10-7 s and 10~‘ s.? ¢; are the precursors of the i-th group of delayed neutrons, 83; 
represents the fraction of the i-th delayed neutron group, and ., is the decay constant 
for the 7-th precursor group. 

Typically, six groups of delayed neutron precursors are taken into account for a 
nuclear reactor simulation. Equation (1) could be programmed in a straightforward way 
but would require more effort using only standard computing elements. This application 
note shows how a specialized computing element can be devised that implements the 
neutron kinetics of a nuclear reactor with reasonable accuracy. 

Figure 1 shows the basic structure of a such a computing element implementing 
the neutron kinetics for either a ?°°U or 73°Pu reactor. The basic circuit consisting of 





1An overview on reactor simulation can be found in [MorRISON 1962]. Details of neutron kinetics and related topics 
are described in depth in [WEINBERG et al. 1958]. An actual implementation of a nuclear plant simulation is shown in 
[FENECH et al. 1973]. 

28 will be ignored in the following — an external neutron source can always be modelled by adding s to the input of 
the network described here. 
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RL and either CL1 or CL2 is that of an integrator with 
M=RL-CL;, 1<3 <2. 


Choosing RL= 100 kQ“ and CL1= 1 nF yields /* = 10~*s. With the same RL a 
CL2= 100 nF results in /* = 10~?s, two suitable and not too unrealistic values for the 
effective neutron lifetime in a nuclear reactor. The /* to be used in a simulation run is 
selected by switch S1. 

The six groups of delayed neutron precursors and resulting neutrons can be modelled 
by series circuits consisting of RU; and CU, for the case of 7°°U and RP; and CP, for 
?39Py respectively. The capacitor values are determined according to* 


Ley 


with 7 being a scaling factor to get reasonable values for the capacitors. In this case 
7 = 10° was chosen. With the capacitances determined this way, the resistors can be 
determined according to 

AiG 

Table 1 shows the \; and (3; for 2°°U and 73°Pu with the corresponding values for 
the resistor-capacitor networks for both cases.* 

The circuit shown in figure 1 can be used directly in an analog computer setup 
although it is advisable to add two electronic switches, one of which connects RL to 
the summing junction of the operational amplifier only when the computer is in operate 
mode, while the other switch discharges the selected capacitor CL1 or CL2 as well as 
the selected network for the delayed neutron groups when the computer is in initial 
condition mode. 





3Cf. [FENECH et al. 1973, p.132]. C*; and R*; denote the capacitor and resistor values for either the 225U or 
239 Pu case. 

4The values A; and §; are according to [TYROR et al. 1970, p. 22]. Actual values tend to differ in the literature a 
bit. 
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Figure 1: Basic circuit for the simulation of neutron kinetics with six groups of delayed neutrons 
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235y 239Dy 

ay Xi Bi C; = 1s Ri [QO] |] A Bi C; = 1s Ri [Q] 

1 |} 0.0127 | 0.000237 | 1.9 ywF 414000 | 0.0129 | 0.0000798 | 619 nF 1252332 

2 || 0.0317 | 0.001385 | 4.37 wF | 72187 0.0311 | 0.000588 1.89 uF 170128 

3 || 0.115 0.001222 | 1.06 wF | 82034 0.134 0.0004536 | 339 nF 220138 

4 || 0.311 0.002645 | 0.85 wF | 37828 0.331 0.0006881 | 208 nF 145248 

5 || 1.4 0.000832 | 59 nF 121065 |) 1.26 0.0002163 | 17 nF 466853 

6 || 3.87 0.000169 | 4.4 nF 587268 || 3.21 0.0000734 | 2.2 nF 1416029 

Table 1: Delayed neutron data for 7°°U and 739Pu reactor 


Due to the rather large values of the series connection of resistors and capacitors, 
the initial condition time should be chosen pretty large to make sure that the capacitors 
are suitably discharged. 


References 


[FENECH et al. 1973] HENRI FENECH, JAY BANE, “A nuclear reactor analog sim- 
ulation for undergraduate nuclear engineering education”, in Simulation, Volume 
20, Issue 4, April 1974, pp. 127-135 


[HUSKEY et al. 1962] HARRY D. HuUSKEY, GRANINO A. KORN, Computer Hand- 
book, McGraw-Hill Book Company, Inc., 1962 


[MORRISON 1962] E. MorRISON, “Nuclear-Reactor Simulation”, in 


[HUSKEY et al. 1962, pp. 5-87-5-93] 


[TYRoR et al. 1970] J. G. Tyror, R. I VAUGHAN, An Introduction to the Neutron 
Kinetics of Nuclear Power Reactors, Pergamon Press, 1970 


[WEINBERG et al. 1958] ALVIN M. WEINBERG, EUGENE P. WIGNER, The Phys- 
ical Theory of Neutron Chain Reactors, The University of Chicago Press, 1958 





Dr. BERND ULMANN, Issue #30, 11-JUN-2021 


